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In the paper a generalization of the (l+l)-dimensional model by Chodos et al [Phys. Rev. 
D61, 045011 (2000)] has been performed to the case of (2+l)-dimensional spacetime. The model 
includes four-fermion interaction both in the fermion-antifermion (or chiral) and fermion-fermion 
(or superconducting) channels. We study temperature T and chemical potential /i induced phase 
transitions in the leading order of large- TV expansion technique, where N is a number of fermion 
fields. It is shown that at sufficiently large values of fi and arbitrary relations between coupling 
constants, superconducting phase appears in the system both at T = and T > 0. In particular, at 
T — and sufficiently weak attractive interaction in the chiral channel, the Cooper pairing occurs 
, for arbitrary couplings in the superconducting channel even at infinitesimal values of p. 
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Last years great attention has been paid to investigation of (2+l)-dimensional quantum field theories (QFT) and, 
in particular, to models with four-fermion interactions of the Gross-Neveu (GN) [l[ type. Partially, this interest 
| is explained by more simple structure of QFT in two-, rather than in three spatial dimensions. As a result, it is 
much easier to investigate qualitatively such real physical phenomena as dynamical symmetry breaking [l|-|8| and 
pi [ color superconductivity @, and to model phase diagrams of real quantum chromodynamics (QCD) [l(| etc. in the 
framework of (2+l)-dimensional models. Another example of this kind is spontaneous chiral symmetry breaking 
induced by external magnetic or chromomagnetic fields. This effect was for the first time studied also in terms of 
(2+l)-dimensional GN model [ll|. Moreover, these theories are very useful in developing new QFT techniques like 
t-H ■ the optimized perturbation theory [l(| [Hf , and so on. 

However, there is yet another more serious motivation for studying (2+l)-dimensional QFT. It is supported by the 
fact that there are many condensed matter systems which, firstly, have a (quasi-)planar structure and, secondly, their 
excitation spectrum is described adequately by relativistic Dirac-like equation rather than by Schrodinger one. Among 
these systems are the high-T c cuprate and iron superconductors i" 1 iij . the one-atom thick layer of carbon atoms, or 
graphene, [hH . [To| etc. Thus, many properties of such condensed matter systems can be explained in the framework 
of various (2+l)-dimensional QFT, including the GN-type models (see, e.g., [i"6U24T | and references therein). 
' In this paper we study phase transitions in a (2+l)-dimensional GN-type model which describes competition 
between two processes: chiral symmetry breaking (excitonic pairing) and superconductivity (Cooper pairing). Clearly, 
the model is suitable for qualitative analysis of superconducting phase transitions in quasi-planar condensed matter 
systems. The structure of our model is a direct generalization of known (l+l)-dimensional model of Chodos et al. 
(2fjl [2fj| . which remarkably mimics the temperature T and chemical potential /i phase diagram of real QCD, to the 
case of (2+l)-dimensional spacetime. Recall that in [25], in order to avoid the prohibition on Cooper pairing as well as 
spontaneous breaking of continuous symmetry in (l+l)-dimensional models (known as the Mermin-Wagner-Coleman 
no-go theorem [27j]), the consideration was performed in the leading order of 1/A^-technique, i.e. in the large- TV limit 
assumption, where N is the number of fermion fields. In this case quantum fluctuations, which would otherwise 
destroy a long-range order corresponding to spontaneous symmetry breaking, are suppressed by 1/N factors. By the 
same reason in (2+l)-dimensional spacetime and in the case of finite values of TV, spontaneous breaking of continuous 
symmetry is allowed only at zero temperature, i.e. it is forbidden at T > 0. Hence, in order to make investigation of 
superconducting phase transitions possible at T > 0, we suppose, as it was done in [25j j . that in the framework of our 
model N — > 00. 

So at T — the results of our paper may be aplied for the description of superconductivity in different TV-layer 
condensed matter systems (N is finite and can even be equal to one), whereas at T > it is better to use the results 
in the description of macroscopic systems composed of a very large number of layers, such as graphite, etc. 

The paper is organized as follows. In Sec. II the GN-type model with four-fermion ineractions in the fermion- 
antifermion (or chiral) and fermion-fermion (or superconducting) channels is presented. Here the unrenormalized 
thermodynamic potential (TDP) of the model is obtained in the leading order of large- N expansion technique. In 
the next Sec. Ill a renormalization group invariant expression for the TDP is obtained whose global minimum point 
provides us with chiral and Cooper pairs condensates. In Sec. IV phase structure of the model is described at T = 
both at fj, = and fi 7^ 0. In particular, it is established in this Section that infinitesimal chemical potential induces the 
superconductivity phenomenon in the case of a rather weak attractive interaction in the fermion-antifermion channel. 
Finally, in Sec. V the (fi, T)-phase diagrams are presented for some representative values of coupling constants. We 
show in this Section that at arbitrary fixed T > superconductivity is induced in the system at sufficiently large 
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values of fi. Some related problems of our consideration are relegated to three Appendices. 



II. 



THE MODEL AND ITS THERMODYNAMIC POTENTIAL 



Our investigation is based on a (2+l)-dimensional GN-type model with massless fermions belonging to a funda- 
mental multiplet of the auxiliary O(N) flavor group. Its Lagrangian describes the interaction both in the scalar 
fermion-antifermion and scalar difermion channels: 
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where fi is the fermion number chemical potential (see also the comments after formula (3)). As noted above, all 
fermion fields tpk (k — 1,...,N) form a fundamental multiplet of O(N) group. Moreover, each field ipk is a four- 
component Dirac spinor (the symbol T denotes the transposition operation). The quantities 7^ {y = 0,1,2) are 
matrices in the 4-dimensional spinor space. Moreover, C = 7 2 is the charge conjugation matrix. The algebra 
of the 7^-matrices as well as their particular representation are given in Appendix [A] Clearly, the Lagrangian L 
is invariant under transformations from the internal auxiliary O(N) group, which is introduced here in order to 
make it possible to perform all the calculations in the framework of the nonperturbative large- AT expansion method. 
Physically more interesting is that the model (1) is invariant under the discrete chiral transformation, tpk — ► 7 5 ^fc (the 
particular realization of the 7 5 -matrix is presented in Appendix [5]) . as well as with respect to the transformations 
from the continuous U(l) fermion number group, "0fc ~~ ^ exp(ia)V>fc (k = 1, N), responsible for the fermion number 
conservation or, equivalently, for the electric charge conservation law in the system under consideration. 

The linearized version of Lagrangian (JTJ that contains auxiliary bosonic fields <r(x), A(x) and A*(x) has the 
following form 
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Clearly, the Lagrangians ([T]) and © are equivalent, as can be seen by using the Euler-Lagrange equations of motion 
for bosonic fields which take the form 



^) = -2^E^' A(x)=-2^^^C^, A*(x) = -2^J2^ 
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One can easily see from ([3]) that the neutral field <r(x) is a real quantity, i.e. (a(x)y = cr(x) (the superscript symbol 
f denotes the Hermitian conjugation), but the (charged) difermion fields A(x) and A*(x) are mutually Hermitian 
conjugated complex quantities, so (A(x))^ = A*(x) and vice versa. Clearly, all the fields © are singlets with respect 
to the auxiliary O(N) group. □ Moreover, with respect to parity transformation P (see also the comment in Appendix 
A), 



P 



ipk(t,x,y) -> 7 5 7Vfc(*, ~x,y), 



1,...,N, 



(4) 



the fields cr(x), A(x) and A*(x) are even quantities, i.e. they are scalars. If the difermion field A(x) has a nonzero 
ground state expectation value, i.e. (A(x)) 7^ 0, the Abelian fermion number £7(1) symmetry of the model is 
spontaneously broken down and the superconducting phase is realized in the model. However, if (<t(x)) 7^ then the 
discrete chiral symmetry of the model is spontaneously broken. 

Let us now study the phase structure of the four- fermion model (1) starting from the equivalent semi-bosonized 
Lagrangian ©. In the leading order of the large- N approximation, the effective action <S e ff (cr, A, A*) of the considered 
model is expressed by means of the path integral over fermion fields 
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1 Note that the A(x) field is a flavor O(N) singlet, since the representations of this group are real. 
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The fermion contribution to the effective action, i.e. the term S c g in ([5]), is given by: 



/ N { N A* A 

JJ[#j][#i] expji / ^ \fik{l"id v + A*7 - <r)i>k - —^C^k ~ ^kC${ d 3 x|. 
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(6) 



The ground state expectation values (a(x)), (A(x)), and (A* (a;)) of the composite bosonic fields are determined by 
the saddle point equations, 



Sa(x) 



= 0, 



6A(x) 



= 0, 



8A*(x) 



= 0. 



(7) 



For simplicity, throughout the paper we suppose that the above mentioned ground state expectation values do not 
depend on space-time coordinates, i.e. 



(a(x))=M, (A(x))=A, (A*(x))^A* 



(8) 



where M, A, A* are constant quantities. In fact, they are coordinates of the global minimum point of the thermo- 
dynamic potential (TDP) D,(M, A, A*). In the leading order of the large- iV expansion it is defined by the following 
expression: 
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d 3 xQ(M,A,A*) = --S eS {a(x),A(x),A*(x)} 



r(x)=M,A(x)=A,A"{x)=A* 



which gives 
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where D = j v id v + ^7° — M. Toproceed, let us first point out that without loss of generality the quantities A, A* 
might be considered as real ones. □ So, in the following we will suppose that A = A* = A, where A is already a real 
quantity. Then, in order to find a convenient expression for the TDP it is necessary to invoke Appendix B, where the 
path integral similar to © is evaluated. So, taking into account in (j9)) the relation (|B7[) we obtain the following 
expression for the zero temperature, T = 0, TDP of the GN model (1): 
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where ) 2 = E 2 + p 2 + A 2 ± 2^M 2 A 2 + p 2 E 2 and E = \J M 2 + |p| 2 . Obviously, the function Q(M, A) is invariant 
under each of the transformations M — > —M, A — > —A and fj, — > —fj,. Hence, without loss of generality, we restrict 
ourselves to the constraints M > 0, A > and fi > and will investigate the global minimum point of the TDP (fTU)) 
just on this region. Using in the expression (110[) a rather general formula 



dp \n (po —A) = yk\A\, 



(11) 



where A is a real quantity, it is possible to reduce it to the following one: 



0(M,A) = n un (M,A) 
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(12) 



The integral term in (|12[) is an ultraviolet divergent one, hence to obtain any information from this expression we 
need to renormalize it. 



2 Otherwise, phases of the complex values A, A* might be eliminated by an appropriate transformation of fermion fields in the path 
integral JgJ. 

3 In Appendix |B] we consider for simplicity the case N = 1, however the procedure is easily generalized to the case with N > 1. 
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III. THE RENORMALIZATION PROCEDURE AT T = 



First of all, let us regularize the zero temperature TDP (|12[) by cutting momenta, i.e. we suppose that < A, 
\p2\ < A in (|T2"j) . As a result we have the following regularized expression (which is finite at finite values of A): 
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Let us use in f)13[) the following asymptotic expansion 
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where \p\ = \/p\ + pf- Then, upon integration there term-by-term, it is possible to find 
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where 0(A°) denotes an expression which is finite in the limit A — > oo. Second, we suppose that the bare coupling 
constants G\ and Gi depends on the cutoff parameter A in such a way that in the limit A — > oo one obtains a finite 
expressions in the square brackets of (jl5[) . Clearly, to fulfil this requirement it is sufficient to require that 
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where are finite and A-independent model parameters with dimensionality of inverse mass. Moreover, since bare 
couplings G\ and G2 do not depend on a normalization point, the same property is also valid for g\.2- Hence, taking 
into account in (|13[) and (|15j) the relations (|16[) and ignoring there an infinite M- and A-independent constant, one 
obtains the following renormalized, i.e. finite, expression for the TDP 
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It should also be mentioned that the TDP (JTTJ) is a renormalization group invariant quantity. 

The fact that it is possible to renormalize the effective potential of the initial model (1) in the leading order of 
the large iV-expansion is the reflection of a more general property of (2+l)-dimensional theories with four-fermion 
interactions. Indeed, it is well known that in the framework of the "naive" perturbation theory (over coupling 
constants) these models are not renormalizable. However, as it was proved in in the framework of nonperturbative 
large iV-technique these models are renormalizable in each order of 1 /iV-expansion. 

In vacuum, i.e. at fj, — 0, the O(A ) term in ([TBI can be calculated explicitly, so we have for the renormalized 
effective potential V(M,A) the expression 
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Now, let us obtain an alternative expression for the renormalized TDP (fT"7|) at /i 7^ 0. For this purpose one can 



rewrite the unrenormalized TDP fl l 
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Since the leading terms of the asymptotic expansion (|14[) do not depend on /i, it is clear that the last integral in 
(1191) is a convergent one. Other terms in (|19[) form the unrenormalized TDP (effective potential) at /1 — which is 



4 Vacuum TDP is usually called effective potential. 
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FIG. 1. The (gi, p2)-phase portrait of the model at n — 
0. The shorthands I, II and III denote the symmetric, 
the chiral symmetry breaking and the superconducting 
phases, respectively. In the phase II (a) = — 1/gi. In the 
phase III (A) = —l/g^. On the line L= {((71,(72) : gi = 
32}, where gi.2 < 0, the TDP minima corresponding to 
the phase II and III are equivalent. 



FIG. 2. The (fi, (72) -phase portrait of the model and crit- 
ical chemical potential fj, C rit(g2) vs gi at arbitrary fixed 
gi < 0. At each point fi = (J, C rit(g2) 7^ there is a first 
order phase transition from the chiral symmetry breaking 
phase II to the superconducting phase III. 



reduced after renormalization procedure to the expression ([18)) . Hence, after renormalization we obtain from (1191) the 
following finite expression (evidently, it coincides with renormalized TDP (1171) ): 



n™(M, a) = v(m, A)- J ^p(eX + s- a - vW+W+Ap - VW+W^W 



(2 

where V(M, A) is presented in (fP5)) . The integral term in ([20)1 can be explicitly calculated. As a result, we have 

(\M 2 fiA 2 / , \3 , 

12TTfl ren (M, A) = + - — + 2 [M + vV + A 2 ) |2M- \/ /1 2 + A 2 



(20) 
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where t± = My/ /j 2 + A 2 ± /i 2 . It is not so evident, but at /i = the expression (|2"Tj) for Q ren (M, A) coincides with 
V(M, A) P 



IV. PHASE STRUCTURE OF THE MODEL AT T = 



As was mentioned above, the coordinates of the global minimum point (Mq,Ao) of the TDP il ren (M, A) define 
the ground state expectation values of auxiliary fields a(x) and A(x). Namely, M = (a(x)) and A = (A(x)). 
The quantities M and A are usually called order parameters, or gaps, because they are responsible for the phase 
structure of the model or, in other words, for the properties of the model ground state (see also the comment after 
(|4|). Moreover, the gap Mq is equal to the dynamical mass of one-fermionic excitations of the ground state. As a 
rule, gaps depend on model parameters as well as on various external factors. In our consideration the gaps Mq and 
Ao are certain functions of the free model parameters <?i and 172 and such external factors as chemical potential \x and 
temperature T. 



A. The case /j, = 



First of all, let us discuss the phase structure of the model (1) in the simplest case when /x = and T = 0. The 
corresponding TDP is given in (I18[) by the function V(M, A). Since the global minimum of this function was already 
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FIG. 3. Superconducting gap Ao = A cr u(g2) vs g-z which 
is generated at the critical point, i.e. at /i = /-i cr it(g2), at 
arbitrary fixed gi < 0. 



FIG. 4. Particle density n = ra C ri*(<?2) vs g2 which is 
generated at the critical point, i.e. at fi = /x C r«(ff2), at 
arbitrary fixed gi < 0. At /i < /x C ri*(<72) the particle 
density n is equal to zero. 



investigated in [28j . although in the framework of another (2+l)-dimensional GN model, we present at once the phase 
structure of the initial model (1) at \i = (see Fig. 1). 

In Fig. 1 the phase portrait of the model is depicted depending on the values of the free model parameters g\ and 
gi. There the plane (31,32) is divided into several areas. In each area one of the phases I, II or III is implemented. 
In the phase I, i.e. at g\ > and gi > 0, the global minimum of the effective potential V(M, A) is arranged at the 
origin. So in this case we have Mq — (u(x)) — and Ao = (A(x)) =0. As a result, in the phase I both discrete chiral 
and continuous electromagnetic U(l) symmetries remain intact and fermions are massless. Due to this reason the 
phase I is called symmetric. In the phase II, which is allowed only for g\ < 0, at the global minimum point (Mo, Ao) 
the relations Mq = — l/gi and Ao = are valid. So in this phase chiral symmetry is spontaneously broken down 
and fermions acquire dynamically the mass M> Finally, in the superconducting phase III, where gi < 0, we have the 
following values for the gaps Mq = and Ao = —l/g%. 

Note also that if g\ = 32 = 9 and, in addition, g < (it is just the line L in Fig. 1), then the effective potential 
(p~8|) has two equivalent global minima. The first one, the point (Mo = — l/g, Ao = 0), corresponds to a phase with 
chiral symmetry breaking. The second one, i.e. the point (Mq — 0, Ao = —l/g), corresponds to superconductivity 

Clearly, if the cutoff parameter A is fixed, then the phase structure of the model can be described in terms of 

bare coupling constants G\ , Gi instead of finite quantities 31 , 32 ■ Indeed, let us first introduce a critical value of the 

2 

couplings, G c = 8A i n (i +x /2) • Then, as it follows from Fig. 1 and (1161) . at G\ < G c and G 2 < G c the symmetric phase 

I of the model is located. If G\ > G C1 Gi < G c (G\ < G c , Gi > G c ), then the chiral symmetry broken phase II 
(the superconducting phase III) is realized. Finally, let us suppose that both G\ > G c and Gi > G c . In this case at 
Gi > G2 (Gi < G2) we have again the chiral symmetry broken phase II (the superconducting phase III). 



B. Consideration of the chemical potential 



In this section we study the influence of the chemical potential /1 > on the phase structure of the model (1) 
(temperature is still vanishing). Numerical and analytical investigations of the TDP (I21[) show that its minimum 
points are of the form (M 7^ 0, A = 0), (M = 0, A ^ 0) or (M = 0, A = 0) only. So to study the properties of the 
global minimum point of the function (|2 1|) it is enough to consider its reductions on the M- and A-axes, where the 
TDP (HO becomes 
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FIG. 5. Superconducting gap Ao and particle density n 
vs fj, at arbitrary fixed gi < and c/2 = 0.5|<?i[. Curves 1 
and 2 are the plots of the dimensionless quantities |(/i|Ao 
and |<?i| 2 n, correspondingly. Here \gi\fic = \gi\licrit(g2 = 
0.5|<7i|) w 0.995 and jgijA. = | 5 i|A orrt (s2 = 0.5| S i|) « 
0.098. 
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FIG. 6. Superconducting gap Ao and particle den- 
sity n vs jj, at arbitrary fixed g± < and $2 = 
— 1.5[<?i|. Curves 1 and 2 are the plots of the dimen- 
sionless quantities |<?i|Ao and |gi| 2 n, respectively. Here 
\gi\fic = |ffi|^cHt(ff2 = -1.5|ffi|) « 0.545 and |ffi|A c = 
|fli]A crtt (52 = -1.5|gi|) « 0.838. 



respectively. Comparing the minima of the functions (|22p and (I23p , it is possible to find the global minimum point of 
the whole TDP (|2ip and its dependence on the model parameters /i, g\ and g 2 , i.e. to determine the phase structure 
of the model. In addition, in the present section we will study the behavior of a particle density n in different phases 
when fi varies, 



dn ren (M, A) 



dfj. 



M=Mo,A=Ao 

coincides with the GMP either of the function u)\ (M) 



(24) 



Since the global minimum point (GMP) of the TDP flU) 
or the function 1^2 (A) (|23|) . it is clear that in the chirally broken phase Ao = and the gap Mo does not depend on 
the parameter g 2 . Correspondingly, in the superconducting phase we have Mq = and the gap Ao does not depend 
on the parameter g%. So, one can use the following expressions for the particle density in the chiral symmetry broken 
II and superconducting III phases: 
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I phase II 



I phase III 
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+ v ^TAg 
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(25) 
(26) 



where 0{x) is the Heaviside step-function. 

The case gi < 0. First of all, let us suppose that g\ is fixed and negative, i.e. g\ < 0. Then it is easy to show that 
for arbitrary value of 52 there exists a critical chemical potential n C rit{92) (see Fig. 2) such that at /i < /i c rtt(<?2) 
the system is in the chiral symmetry breaking phase II (if fi C rit(g2) > 0), and it is in the superconducting phase III 
at fi > (x C rit(g2)- In other words, if /1 < fi C rit(92) 7^ 0, then the global minimum of the TDP (|21j) lies at the point 
(Mo = — l/<7i,Ao = 0) which does not depend on fi in the interval < fi < fi C rit(g2)- However, at /x = fJ, C rit{g2) 
it jumps to the point (Mo = 0, Ao = A cr i t (g2)), where A cr it(<?2) vs 32 is depicted in Fig. 3. Hence, at the critical 
point \i = /i cr it(<?2) a first order phase transition occurs and a superconducting gap Ao = A cr it(<?2) is dynamically 
generated. It turns out that A vs n is an increasing function in the interval \x > fJ. cr it{92)- In particular, the behavior 
A vs fj, is presented in Fig. 5 (at g 2 — 0.5|<7i|), Fig. 6 (at g 2 = — 1.5|gi|) and Fig. 7 (at g 2 = —0.5\gi\) as the curve 1. 

Moreover, it is clear from Fig. 2 that at < fJ, C rit{92), i.e. in the phase II, the particle density is equal to zero. To 
explain this circumstance, recall that in the phase II the gap Mo is equal to l/|gi|. So, for all g 2 -vahies the relation 
H>crit{g2) < Mo is valid (see Fig. 2). As a result, throughout the phase II, where fi < fi C rit{g2), we have fi < Mq and 
hence, as it follows from the relation (|25p . the zero particle density, n = 0. However, when fi reaches its critical value, 
M = Merit (32)) the nonzero particle density n cr .;t(<?2) is generated dynamically in the system (see Fig. 4). Further 



5 All the Figs 2-10 are drawn in terms of dimensionless quantities which arc obtained after multiplication of appropriate powers of |gi| 
with corresponding dimensional quantities. For example, there instead of fi, Ao, gi we use their dimensionless analogies |gi|^t, |gi|Ao, 
52 / 131 1 • Instead of particle density n the dimensionless quantity g\n is depicted there etc. 
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FIG. 7. Superconducting gap Ao and particle density n vs 
fi at arbitrary fixed gi (both at gi < and gi > 0) as well 
as at <?2 = — 0.5 1 <7i | . Curves 1 and 2 are the plots of the 
dimensionless quantities |<?i|Ao and 6|<;i| 2 n, respectively. 



FIG. 8. Superconducting gap Ao and particle density n 
vs /j, at arbitrary fixed gi > as well as at gi — 0.5<?i. 
Curves 1 and 2 are the plots of the dimensionless quanti- 
ties |<7i|Ao and |(?i| 2 n, respectively. 



growth of the chemical potential is accompanied by increase of the particle density n vs //. (Evidently, in this case 
the particle density must be calculated with the help of the expression (J2U).) For example, in Figs 5-7 at the same 
representative relations between g\ and 172 the particle density n vs /1 is depicted as a monotonically increasing curve 
2. 

Finally recall that at /1 = the two phases, II and III, have equivalent minima of the TDP only at negative values 
of gi — gi (it is the line L in Fig. 1). It turns out that for arbitrary fixed g\ < and at growing chemical potential, 
this property of the TDP is also allowed but in a much more extensive c/2-region. Indeed, as our analysis shows in this 
case, if 172 > or 52 < gi then at /j, = /i cr it(ff2) (see Fig. 2) the TDP has two equivalent minima, corresponding to 
these phases. As a result, for these values of 31 and gi there is a coexistence of chirally broken and superconducting 
phases at fi = (J, C rit{92)- In this case, when viewed from the side we have the following picture of phase transitions in 
the system. At rather small values of /i the ground state of the system is an empty space (particle density is zero). 
If fermions are created in this state, they have a mass equal to Mq = —1/gi, i.e. the ground state corresponds to a 
chirally broken phase II. Then, if /i reaches the critical value [i = \i C Ht (92), bubbles of a new phase III appear in the 
empty space. Inside each bubble the particle density n is nonzero and equal to n cr it(gi) (see Fig. 4). 

The case g\ > 0. Now the model phase structure consideration for a positive gi-values is in order. Recall, in 
this case we have a rather weak attractive interaction in the chiral channel, i.e. G\ < G c . Evidently, if in addition 
gi < 0, then in this case the superconducting phase is realized for arbitrary values of fi > 0. The behavior of the gap 
Ao and particle density n vs [i in this branch of the superconducting phase is given in Fig. 7 in the particular case 
gi = — 0.5<7i for g\ > 0. Moreover, as it is clear from Fig. 7, the same behavior for Ao and nvs/i remains valid for 
the case 172 = — 0.5|gi| and negative values of g\. To explain this fact, it is necessary to take into account the remark 
made after formula (|24[) that the superconducting gap does not depend on the coupling 171 but only on gi one. So, it 
is no wonder that the plots of Ao and n are not changed when the parameter gi changes the sign. 

Recall, if both g\ > and gi > 0, then we have at /1 = the phase I without any symmetry breaking, where the 
gaps Ao and Mq vanishes (see Fig. 1). However, our analysis shows that at arbitrary small nonzero /i the global 
minimum point of the TDP (f2Tj) moves from the point (Mo = 0, Ao = 0) to the following one (Mo = 0, Ao 7^ 0). 
Hence, at positive values of g\ and gi a continuous second order phase transition occurs from symmetric phase I to 
superconducting one III when chemical potential acquires an arbitrary small nonzero value. The typical behavior of 
the gap Ao and particle density n vs fi in this superconductivity region is depicted in Fig. 8. Comparing Figs. 7 and 
8, we see that at the same value of [i the gap Ao and particle density n are much greater in the case 171 > 0, gi < 0, 
than in the case g\ > 0, gi > 0. To support this statement we draw in Figs. 9 and 10 the plots of the gap Ao and 
particle density n vs gi in two different regions gi < and gi > 0, respectively, at the particular value of the chemical 
potential, fi — 0.5/ 171. 

We see that at gi > 0, i.e. at G\ < G c , the chiral symmetry breaking is absent but the Cooper pairing phase occurs 
at any /1 > 0. To explain this different behavior, one can use the following very naive physical arguments. Since at 
fi > we have a nonzero particle density (see, e.g., in Fig. 8), there is a Fermi sea of particles with energies less or 
equal to fi (Fermi surface). Evidently, in this case there is no energy cost for creating a pair of particles with opposite 
momenta just over the Fermi surface. Then, due to an arbitrary weak attraction between these particles (Gi > 0), 
the Cooper pair is formed and U(X) symmetry is spontaneously broken, as a result of Bose-Einstein condensation of 
Cooper pairs. Note, since in the energy spectrum of fermions the gap A / appears (see in Fig. 8), rather small 
external forces are not able to destroy the superconducting condensate and it is a stable one. 
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FIG. 9. Superconducting gap Ao and particle density n vs 
Q2 < at arbitrary fixed gi > and fj, — 0.5/<?i. Curves 
1 and 2 are the plots of the dimensionless quantities giAo 
and gin, respectively. 



11.14 - 




FIG. 10. Superconducting gap Ao and particle density 
n vs 92 > at arbitrary fixed gi > and /j, — 0.5 jg\. 
Curves 1 and 2 are the plots of the dimensionless quanti- 
ties <?iAq and gin, respectively. 



Concerning the chiral symmetry breaking in this case, it is clear that a particle and a hole with opposite momenta 
can also be created without any energy cost in the system. Moreover, there is also an attraction between a particle 
and a hole. However, since the nonzero gap M does not appear in the energy spectrum at sufficiently small G\ < G c , 
the particle-hole pairing in this case is rather a weakly bounded resonance, which, unlike a stable pair, could be easily 
destroyed by an arbitrary small external influence. So, no stable Bose-Einstein condensate of these pairs does appear 
and chiral symmetry remains intact. (For a more detailed discussion on possible types of pairing in dense (quark) 
fermionic matter see, e.g., in |29j.) 

In summary, we can say that at T = chemical potential induces superconductivity in the model for arbitrary 
relations between coupling constants g±.2 (or, equivalently, Gi^)- 



V. FINITE TEMPERATURE 



Now let us study the influence of both temperature T and chemical potential /i on the phase structure of the model. 
It is well known (see, e.g., in (30j ) that in d space dimensions (in our case, evidently, d = 2) the transition probability 
from one degenerated minimum of the TDP to another is proportional to exp(— N j3L d ~ 2 ), where L is the linear size 
of the system and f3 is the inverse temperature, f3 = 1/T. It follows from this expression that at d = 2 the transition 
probability is zero even at finite N if T = 0. This leads to the fact that a continuous symmetry can be spontaneously 
broken in any planar systems at T = 0. (Hence, our consideration of superconducting phase transitions performed 
at T — in the previous section is valid for arbitrary values of iV.) However, if T ^ 0, then transition probability 
in the above expression does not vanish at finite N . This circumstance ensures vanishing of the order parameter 
and, as a result, might lead to a prohibition for spontaneous symmetry breaking in d = 2 spatial dimensions at finite 
N and T 0. However, if N — > oo the transition probability vanishes and the spontaneous symmetry breaking is 
allowed. Just this assumption, i.e the same as in [25l [26j . is used in the following consideration, where we study the 
temperature dependent superconducting phase transitions in the leading order of large- N expansion technique. 

In this case, in order to get the corresponding (unrenormalized) thermodynamic potential f2 T (M, A) one can simply 
start from the expression for the TDP at zero temperature (|10[) and perform the following standard replacements: 

oo 

J^(...)-HT £ (••■), po->pon = iw n = i7rT(2n+l), n = 0,±1,±2, (27) 

~ n— — oo 

i.e. the po-ffitegration should be replaced by the summation over Matsubara frequencies uj n . Summing over Matsubara 
frequencies in the obtained expression (the corresponding technique is presented, e.g., in [3l|), one can find for the 
TDP 

^ a > - € + S - £ I? k + « - 2T £ w? " d 1 + ^ i 1 + '"^) ■ ,28) 

where /3 = 1/T and are given in (|10[) . Clearly, only the first integral in this expression (which is the same as in 
the zero temperature case) is responsible for ultraviolet divergency of the whole TDP (|2"8|) . So, regularizing the TDP 
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FIG. 11. (fi, T)-phase diagram of the model at gi = 
— 0.5|<7i| and arbitrary fixed g\ both at gi < and gi > 0. 



FIG. 12. (fi, T)-phase diagram of the model at arbitrary 
fixed gi > and at 32 = 0.5<?i. 



in the way as it was done in (fT5j) for zero temperature TDP and then replacing G± t 2 —> Gi,2(A) (see formula 
([in])), we can obtain in the limit A->ooa finite expression denoted as f2£ e ™(M, A), 



CL r T en (M, A) = Q ren (M, A) — 2T 



d 2 p 
-co (2^F 



In [1 



-pei 



[1 



(29) 



where fl ren (M, A) is the zero temperature TDP ([2T]) . Numerical investigations show that all possible local minima 
of the TDP f^ erl (M, A) are located in the lines M = or A = 0. So it is sufficient to deal with corresponding 
restrictions of the TDP on these lines, i.e. with the following functions, 



Fi(M) = n r T en (M, A) 



F a (A) = n r T en {M,A) 



A=0 



M 



(M) - 2T J°° ^ In ( [1 + e -K E+ ri] [l + e^^'l] ' 

(27r)2 ^([l + e-^][l + e -^]), (30) 



A/ 3 
"3V 



M=0 



2tt5i 
w 2 (A)-2T 



/■°° d 2 D 

00 j2 



-00 (2tt) 



1 



(31) 



where F = ^/Ip] 2 + M 2 , (F^) 2 = ± ^i) 2 + A 2 , and the functions wi(M), w 2 (A) are presented in (22) and ([25]). 
respectively. The gaps Mo and Aq are the solutions of the following stationary (gap) equations, 



where 



dF l (M) _M dF 2 (A) _ A 

flu = Jl(-M) = 0, a a =— 72(A =0, 

am tt aA n 



fi(M) = j + M + Tin { [1 + e -^ /+ ">] [1 + e -/3(M-^)] j ^ 
/ 2 (A) = 1 + Vm 2 + A 2 + 2Fln f 1 + e -Py/^^\ - M f tanh f ^ 



x/q 2 + A 2 



(32) 

(33) 
(34) 



respectively (for details, see Appendix [C]) . On the basis of these gap equations we will study the phase structure of 
the model at T > 0. 

The case g\ > 0. First of all let us consider the phase portrait of the model at g\ > 0. It straightforwardly follows 
from (|32p and (155]) that the gap Mq is always zero at g\ > (it is a nonzero quantity only at gi < 0). However, the 
gap A is positive both at g 2 < and g 2 > if temperature is sufficiently small, i.e. when /2(0) < 0. So, at gi > 
and for arbitrary values of fj, the superconducting phase III is arranged at sufficiently small values of temperature 
T < T c (fj). At T > T c (/i) the gap equations (1521) supply the Ao = and Mq = gap values, i.e. the symmetric phase. 
The second order phase transition temperature T c (jj) is the solution of the equation / 2 (0) = 0, 



/ 2 (0) 



1 

.92 



/i + 2Fln(l + e 



-Pv\ 



(35) 
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FIG. 13. (/i, T)-phase diagram of the model at 32 = 
0.5|<7i| and arbitrary fixed gi < 0. All the curves are the 
lines of second order phase transitions except the bound- 
ary between the superconducting and chiral symmetry 
breaking phases, where a first order phase transition is 
realized. The coordinates of the tricritical point A are 
the following ones, |<;i|/m ~ 0.999 and |<?i|Ta « 0.056. 
Moreover, \gi\fi c « 0.995 and |fli|T c = 1/(2 In 2) w 0.721. 




FIG. 14. (/i, T)-phase diagram of the model at gi = 
— 1.5 1 gi I and arbitrary fixed gi < 0. All the curves are the 
lines of second order phase transitions except the bound- 
ary between the superconducting and chiral symmetry 
breaking phases, where a first order phase transition is 
realized. The coordinates of the tricritical point A arc 
the following ones, ~ 0.645 and |<?i|Tk ~ 0.602. 

Moreover, |gi|/x c w 0.545 and |#i|T c = 1/(2 In 2) « 0.721. 



Hence, in the (/x,T)-plane the curve T = T c (/x) is the boundary between symmetric and superconducting phases. 
Numerical investigation of the equation ([55)1 produces at g 2 = ±0.5gi the phase portraits of the model presented in 
Figs 11 and 12. 

If g 2 < 0, then it follows from (|35|) that T c (0) = — l/(2g 2 ln2) . Moreover, in this case the critical temperature can 
be given as a series over the small parameter [i, 

T c (ji) = T c (0) - M 2 <?2/16 + o(p 2 g 2 ). (36) 

Comparing this expansion at g 2 = — 0.5gi with T c (/x) of Fig. 11, we see that (J36J) supplies a rather good approximation 
for the critical temperature only in the interval < \xg\ < 0.2. 

Now let us try to present some analytical approximation for the T c (/i) at 52 > (gi is still fixed and positive). For 
this purpose note first of all that for all points of the critical curve T = T c (fi) of Fig. 12 the relation fi/T = fi(3 >> 1 
is valid. Then, it is convenient to present the equation (|35l) in the following equivalent form: 

^~ + f + In (1 + e-/*) -%U+ r 2 ^ + C a - f° [tanh, - 1]^U 0, (37) 

where 

r 1 dz f°° dz 

Ci= tanhz— « 0.910, C 2 = / [tanh z - 11 — ps -0.091. (38) 
Jo z Ji z 

The third term in ([37]) as well as the last integral in the braces of (|37|) can be neglected in comparison with other 
terms. The obtained equation can be easily solved with respect to T. As a result we have 

T c (jjl) ps I exp [Cx + C 2 - 1 - 1/( M <? 2 )] . (39) 

Note that at g 2 = 0.5<7i the plot of the expression (|39[) coincides with great accuracy with the critical temperature of 
Fig. 12 in the whole interval < figi < 2. 

The case gi < 0. In this case we present three (/i, T)-phase portraits of the model for qualitatively distinct 
relations between g\ and g 2 - The first one for g 2 — — 0.5|<?i| (which is in Fig. 11) was already described above because 
it is the same as in the case g\ > 0, g 2 = —O.bgi. The other two phase portraits are represented in Figs. 13 and 14 
for g 2 = 0.5|<7i| and g 2 = — 1.5|gi|, respectively. There the points (p,T) of the boundary between the symmetric and 
chiral symmetry breaking (or superconducting) phases are given implicitly by the equation /i(0) = (or f 2 (0) = 0), 
where the functions fi(M) and /2(A) are defined in (|33p and (|3"4"1) . respectively. On these boundaries, the second 
order phase transitions occur. In contrast, the boundary between chiral symmetry breaking and superconducting 
phases is the curve of the first order phase transitions. So at the points (/i, T) of this boundary the two phases may 
coexist. 
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Analyzing the cited above (//, T)-phase diagrams of Figs. 11-14, we see that for each arbitrary fixed value T of the 
temperature (and for all relations between coupling constants) there exist a definite value fpp of the chemical potential 
such that for all \i > [It the superconducting phase is realized in the system. This property is inherent only to a 
(2+l)-dimensional model (1) and it is absent in the two-dimensional analogue [25j . 



VI. SUMMARY AND CONCLUSIONS 



In this paper we study the competition between chiral and superconducting condensations in the framework of the 
(2+l)-dimensional GN-type model (1) which is a direct generalization of the two-dimensional analogue by Chodos 
et al. 25]. So, the initial four-fermion model (1) describes interactions both in the fermion-antifermion (or chiral) 
and superconducting difermion (or Cooper pairing) channels with couplings G\ and G2, respectively. Moreover, it is 
chirally and U(l) invariant one (the last group corresponds to conservation of the fermion number or electric charge 
of the system). To avoid the ban on the spontaneous breaking of continuous symmetry in (2+l)-dimensional field 
theories at T > 0, we consider, as it was done in [251 ], the phase structure of our model in the leading order of the 
large- iV technique, i.e. in the limit N — > 00, where N is a number of fermion fields. 

The case T = 0, fi = 0. First of all we have investigated the thermodynamic potential of the model at T = 0, 
fj, = 0. In this case the phase portrait is presented in Fig. 1 in terms of the renormalization group invariant finite 
coupling constants g% and gi. Each point (51,52) of this diagram corresponds to a definite phase. For example, at 
51,2 > 0, i.e. at sufficiently small values of the bare coupling constants (see the comment at the end of Section 
IIV A[) . neither chiral nor U(l) symmetries are violated and the system is in the symmetric phase, etc. 

The case T = 0, /i 7^ 0. In this case we select two qualitatively different situations, g\ < and g\ > 0. If g\ < 
and fixed, then in Fig. 2 we draw the (32, /x)-phase diagram of the model. It means that at 52 > or at gi < g\ 
the phase II with zero particle density is realized at sufficiently low values of /i. In this case the ground state of the 
system is an empty space. Then at some critical value /1 = /UcriiGfe) bubbles of the new phase III with particle density 
n C ru{92) (see Fig. 4) can appear in the space, and for all /1 > (J, C rit(92) the whole space is filled with superconducting 
phase, in which particle density n is not zero, n > n cr it{g2)- If <7i > 0, then the system is in the superconducting phase 
even at arbitrary small values of fi. Hence, at T = and at growing chemical potential the system is transformed 
into a superconducting state. 

The case T > 0, /i 7^ 0. Phase portraits of the model are presented in this case in Figs 11-14. It is clear from the 
figures that at fixed /1 and increasing temperature the symmetric phase is restored. However, at arbitrary fixed T, 
growth of the chemical potential leads to appearing of superconductivity in the system at arbitrary relations between 
coupling constants g\ and gi. 

The fact that chemical potential induces superconductivity phenomenon is the main result of our paper. Note that 
in general this property of the (2+l)-dimensional GN-type model (1) is not valid in the case of the two-dimensional 
model [H. 

We hope that our investigations can shed new light on the superconducting phenomena in condensed matter systems 
with planar structures. 



Appendix A: Algebra of the 7-matrices in the case of SO(2,l) group 

The two-dimensional irreducible representation of the 3-dimensional Lorentz group SO(2,l) is realized by the fol- 
lowing 2x2 7-matrices: 

?="=(i!).*'^=(:;).?=-i'i). < ai > 

acting on two-component Dirac spinors. 
They have the properties: 

Tr(rr) = 1gT\ lT,r] = -2*e^ a ; ^ = -i £ ^j a +g^, (A2) 

where g^ = g^ v = diag(l, — 1, —1), 7 Q = 5a/37' 3 , £ 012 = 1. There is also the relation: 

Tr{Yl v T) = -2ie^ a . (A3) 

Note that the definition of chiral symmetry is slightly unusual in three dimensions (spin is here a pseudoscalar rather 
than a (axial) vector). The formal reason is simply that there exists no other 2x2 matrix anticommuting with the 
Dirac matrices j" which would allow the introduction of a 7 5 -matrix in the irreducible representation. The important 
concept of 'chiral' symmetries and their breakdown by mass terms can nevertheless be realized also in the framework 
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of (2+l)-dimensional quantum field theories by considering a four-component reducible representation for Dirac fields. 
In this case the Dirac spinors ip have the following form: 

= ( tfl) , (A4) 



with ipi,i/>2 being two-component spinors. In the reducible four-dimensional spinor representation one deals with 
(4x4) 7-matrices: 7^ = diag(^^, — 7 /1 ), where 7^ are given in (|Aip . One can easily show, that (/x, v = 0, 1, 2): 

Tr(7^)=4g^; 7^ = ^+3^; 

^ - ^,7*1 = diag(-ie^ a ,-ie^ a %). (A5) 

In addition to the Dirac matrices 7 M (fj, = 0, 1, 2) there exist two other matrices 7 3 , 7 s which anticommute with all 
7 M (/J, = 0,1,2) and with themselves 

^ 3 =(?;0' 75=7 ° 7V73=i (?;"o / )' (A6) 

with I being the unit 2x2 matrix. Finally note that in terms of two-component spinors tpi , ip2 the parity transformation 
P, defined in the space of four-component spinors by the relation Q, looks like 

P : 4>i(t,x, y) -> rffoit, -x,y); ^(t, x, y) -> «7 1 V'i(t, -x, y). (A7) 

Such a definition of the space parity transformation is commonly used in (2+l)-dimensional theories with four- 
component representation for Dirac spinors (see, e.g., in [32l|). 



Appendix B: The path integration over anticommutating fields 



/ = / [dq][dq] exp (i / d 3 x[qDq - -(q T Cq) - -(qCq T )] ) , 



(Bl) 



Let us calculate the following path integral over anticommutating four-component Dirac spinor fields q(x), q{x): 

A , rr_ , A 

2" 

where we use the notations of Section HTl and, in particular, the operator D is given in (|9]). Note in addition, the 
integral I is equal to the argument of the In- function in the formula (|9|) in the particular case N = 1. Recall, there 
are general Gaussian path integrals |33j : 

' (' \<i ■- 11 1 11 i = [(let ( ill' " •-<<; ! - > 

i 



[dq] exp (i 
[dq] exp 



d 3 3 



-q T Aq + 



(det(A)) 1/2 exp(- ' 



d 3 a 



d 3 a 



l -qAf + m T )= (det(A)) 1 / 2 exp ( - I / d 3 x fjA' 1 ^ 



1 — 1 -T 



(B2) 
(B3) 



where A is an antisymmetric operator in coordinate and spinor spaces, and T)(x), fj(x) are anticommutating spinor 
sources which also anticommutate with q and q. First, let us integrate in (|B1I) over g-ficlds with the help of the 
relation (|B2I) supposing there that A = AC, qD — rj T , i.e. r\ — D T q T . Then 



1= (dct(AC)) 1/2 J[dq]exp( - % - J d 3 xq[AC + D(AC)- 1 D?]? 



(B4) 

Second, the integration over g-fields in (|B4|) can be easily performed with the help of the formula (|B3j) . where one 
should put A = AC + D(ACy 1 D T and fj — 0. As a result, we have 



I = (dct(AC)) 1/2 (Aet[AC + DiACy^T]) 112 = (det[A 2 C 2 + DC' 1 D T C]) 
Taking into account the relations {d v ) T = —d u and C~ 1 ("f") T C = — 7^ (v = 0, 1, 2), we obtain from (|B5 



,V2 



I = (det[-A 2 + D + D- 



,1/2 



(det B) 



1/2 



where D± = 7 l/ i<9„ — M ± //7 . Using the general relation det-B = exp(Trln_B), we get from 



In/ = ^Trln (B) = £ J -^L H X l{p) ) J <f x . 



(B5) 
(B6) 

(B7) 



(A more detailed consideration of operator traces is presented in Appendix A of the paper |34[.) In this formula 
symbol Tr means the trace of an operator both in the coordinate and internal spaces. Moreover, Xi(p) (i = 1,2) in 
(|B7|) are two twice degenerated eigenvalues of the 4x4 Fourier transformation matrix B(p) of the operator B, i.e. 



Ai, 2 (p) = M 2 



2 2 2,2 
Pi - P2 ~ A 1 + Po 



A 2 ± 2J-M 2 p\ - M 2 p\ + M 2 pl + p?p\ + p?p\. 



(B8) 
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Appendix C: Gap equations 

The equation for the gap Mo, i.e. the first one of equations (|32|) . is obtained, e.g., in [4j, where a phase structure 
of the initial model (1) was consided in the particular case of G 2 — 0. 

To obtain a gap equation for the superconducting gap A , 0F 2 (A)/dA = 0, let us first transform the original 
expression (l3lT) for the TDP F 2 (A) using polar coordinates in the integral in (l3i~j) . Integrating in the obtained 
expression over a polar angle, we have 

TpOO rp 



(A) = wa(A) - - j pdpln (l + e -/Vb+M) 2 +A 2 ) _ £ j pdpln ^ + e -^( P -^+A^ (C1) 

It is very convenient to change integration variables in (|C1|) (we use q = p + [j, for the first integral and q = p — \i for 
the second one) and, after some manipulations, to get an equivalent expression, 

F 2 (A) = wa(A) - ^ ^ Qdq In (l + e^V^+A 5 ) _ 2 J± £ dq m (i + e -fi sf¥T&^ . ( C 2) 
Starting from (|C2[) and taking into account the expression (j2"5)l for ^(A), we have the following gap equation: 



OA 7TI72 7T 7T \ / 

2qdg 2A/i /"^ 



The first integral in (|C3[) is a rather simple one, i.e. 



0. (C3) 



2qdq 



^/q 2 + A 2 (l + e ^i 2 +^- 



In ( 1 + e -^V^ 2 +A 2 ) . (C4) 



In contrast, let us present the third term in (|C3[) in the integral form, i.e. 



.a ln wjTAn = _^ ^ ^ 



which then can be combined with the last integral of (|C3|) . As a result we obtain for the superconducting gap Ao the 
second of equations (|32|) . 
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